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Theoretical description of nuclear fission remains one of the major challenges of quantum many- 
body dynamics. The slow, mostly adiabatic motion through the fission barrier is followed by a 
fast, non-adiabatic descent of the potential between the fragments. The latter stage is essentially 
unexplored. However, it is crucial as it generates most of the excitation energy in the fragments. The 
superfluid dynamics in the latter stage of fission is obtained with the time-dependent Hartree-Fock 
theory including BCS dynamical pairing correlations. The fission modes of the ^^®Fm nucleus are 
studied. The resulting fission fragment characteristics show a good agreement with experimental 
data. Quantum shell effects are shown to play a crucial role in the dynamics and formation of the 
fragments. The importance of quantum fluctuations beyond the independent particle/quasi-particle 
picture is underlined and qualitatively studied. 
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Since its discovery in 1939 HU], nuclear fission has 
been a pillar of nuclear physics. It is indeed the pro¬ 
cess revealing most clearly the complexity of low-energy 
nuclear dynamics. Thus, it provides an ideal test for 
the modelling of nuclear systems with quantum many- 
body theories. It has also found important applications 
in fundamental science. For instance, it is the best way 
to produce beams of exotic rare isotopes in present and 
future accelerators. In addition, understanding fission 
dynamics is crucial for the production of superheavy el¬ 
ements, an important motivation for the construction of 
exotic beam accelerators. Fission is also present in some 
astrophysical processes. In fact, the natural abundance 
of elements heavier than iron is believed to be largely 
influenced by the fission of very neutron-rich heavy nu¬ 
clei formed in supernovae and neutron star mergers [3]. 

Moreover, nuclear fission is one of the greatest sources of 
energy available on Earth. Safely extracting this energy 
to produce electricity while preserving our environment 
has been, and still is, one of the greatest technological 
challenges of humanity. 

The theoretical description of the nuclear fission phe¬ 
nomenon remains a profound problem in fundamental 
science. In order for fission to occur, nuclei have to over¬ 
come the fission barrier, which involve slow, dissipative 
motion. This process has been widely modelled using 
the adiabatic approximation [4]. This approximation is 
justified near the fission barrier by the fact that the evo¬ 
lution of the collective coordinate is slow enough to allow 
the internal degrees of freedom to be equilibrated. How¬ 
ever, after the barrier is passed, the fission fragments un¬ 


dergo an accelerating descent to scission, which involves 
non-adiabatic effects, and so cannot be described with 
standard adiabatic models. Thus, the adiabatic approx¬ 
imation breaks down in the latter stage of fission, where 
the evolution is faster. In particular, the dynamics near 
scission, where the fragments separate, is clearly non- 
adiabatic H Ej. These effects are crucial to properly 
describe properties of the fragments such as their mass, 
charge, and their excitation energy. In particular, the 
latter determines the number of emitted neutrons and is 
thus one of the most important properties for the simu¬ 
lation and safety of future nuclear reactors. 

The complexity of fission dynamics and the high num¬ 
ber of degrees of freedom to be included motivate the use 
of microscopic approaches, where the quantum behaviour 
of each and every nucleon in the whole system is followed 
in time. Microscopic approaches have recently consider¬ 
ably improved our understanding of the fission process 
IMS]. Time dependence is also a key to addressing the 
latter stage of the fission process [20]. Non-adiabatic ef¬ 
fects in the latter stage of fission have then recently been 
investigated in fission of ^^^Fm using a time-dependent 
mean-field approach [18]. As an example of new out¬ 
comes, it was shown that more than half of the final 
excitation energy of the fragments is acquired during the 
last zeptosecond before scission and that it is at least 
partly stored in low-energy collective vibrational states 
of the fragments. However, the calculations in Ref. m 
are based on an independent particle approximation and 
pairing correlations responsible for a superfluid phase in 
nuclear systems were neglected, thus limiting the range 
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of possible applications to very few non-superfluid sys¬ 
tems. It is then crucial to incorporate time-dependent 
pairing correlations in order to investigate fission dynam¬ 
ics across the nuclear chart. These correlations have been 
recently included in realistic time-dependent mean-field 
calculations [2TH25] which we extend to the study of fis¬ 
sion. 

The purpose of this letter is to present a microscopic 
method which incorporates both superfluid dynamics and 
non-adiabatic effects in the latter stage of the fission pro¬ 
cess. Fission in the ^^^Fm nucleus is considered as an 
example of application. This nucleus is known exper¬ 
imentally to exhibit a bimodal fission [26] and consti¬ 
tute an ideal benchmark for theoretical studies of fission 

OEunmiiiiH]. 

In the present approach, it is assumed that the fission 
process is divided in two steps. In a first step, the slow 
evolution near the fission barrier is treated in a standard 
way using the adiabatic approximation. The constrained 
Hartree-Fock (HF) equations with pairing correlations 
are solved at the BCS level (CHF+BCS). Several con¬ 
straints are considered in order to find different valleys 
in the potential energy surface. In a second step, the 
non-adiabatic descent of the potential towards scission is 
determined using the time-dependent HF equations with 
dynamical pairing correlations (TDHF+BCS). The prop¬ 
erties of the fragments, in particular their mass, charge, 
and kinetic energy, are then computed after scission and 
compared with experimental data from [26] . 
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FIG. 1: (Color online) Potential energy in the three valleys: 
symmetric compact fragment (scf) (blue dashed line), sym¬ 
metric elongated fragment (sef) (solid red line) and asymmet¬ 
ric elongated fragment (aef) (green doted line). The arrows 
correspond to the starting configuration of the dynamical cal¬ 
culations for each mode. Isodensities at half the saturation 
density po/2 = 0.08 fm“^ are also plotted at these initial 
configurations and for the ^^®Fm ground-state (left). 

The mean-field is obtained with the Slydd [29] Skyrme 
energy density functional and a constant-G interaction in 
the pairing channel. The pairing strength is chosen to be 
Gn = 24/Ar MeV and Gp = 19/Z MeV. The CHF+BCS 
and TDHF+BCS calculations are obtained with modified 


versions of the ev8 m and tdhfSd [29] codes, respec¬ 
tively, assuming only one plane of symmetry. All calcula¬ 
tions are performed on a Cartesian grid of 88x19.2x19.2 
fm^ with a mesh size 0.8 fm. The time evolution is ob¬ 
tained with a time step 1.5 x 10“^^ s. 

The CHF+BCS solutions with constraints on 
quadrupole Q 2 o^ octupole Qso, and hexadecapole Q^o 
moments along the fission axis let appear three valleys. 
Two symmetric valleys (i.e., with a total Qsq = 0) lead 
to symmetric compact fragments (scf) and to symmetric 
elongated fragments (sef), respectively. In the scf valley, 
the final fragments are almost spherical while they 
exhibit a strong prolate shape in the sef valley. A third 
valley with Qso 7^ 0 leads to asymmetric elongated 
fragments (aef) with different masses and charges. 
Similar adiabatic valleys were obtained by other groups 
for the same nucleus m 0 ng. The potential energy 
along these three valleys is shown in Fig. This 

first, adiabatic, stage of fission is crucial in determining 
the outcome of the reaction. We see that the scf and 
aef valleys have similar energies up to relatively large 
deformations {Q 20 ^ 270 fm^) in the descent to scission. 
The sef valley, however, is found at higher energy. 

Let us now investigate the second stage of the fis¬ 
sion process, associated with the non-adiabatic descent 
of the potential towards scission. TDHF+BCS calcu¬ 
lations have been performed with initial configurations 
along theses valleys indicated by arrows in Fig. More 
compact configurations belong to the adiabatic phase as 
single particle states cross the Fermi level and induce 
a two-body dissipation due to the Landau-Zener effect 
[n in]. The density evolutions in the non-adiabatic 
phase are represented at various times for each mode in 
Fig. [21 The asymmetric mode is likely to be responsible 
for the tail of the experimental fragment mass distribu¬ 
tion shown in the inset of Fig. [^ It is also interesting 
to note that the three evolutions require different times 
to reach scission. These times are ^ 2 zs, ^ 5.4 zs, and 
^ 3.2 zs for the scf, aef, and sef modes, respectively. 
These variations are likely to be due to a combination of 
two factors: different potential slopes and different one- 
body viscosity which is expected to depend on shell ef¬ 
fects. It is worth mentioning that the dynamical pairing 
effects significantly affect the fission dynamics. Indeed, 
freezing pairing correlations by keeping occupation num¬ 
bers constant in time, the so-called Frozen Occupation 
Approximation (FOA) [25] [32] , inhibits strongly the fis¬ 
sion process. 

The final total kinetic energy (TKE) of the fragments 
is another important observable which can be used to 
distinguish between the fission modes. In purely adia¬ 
batic approaches, the TKE is usually estimated from the 
scission configuration which is identified on the potential 
energy surface, based on some criteria [6] |8j [19] . The ad¬ 
vantage of using the TDHE-FBCS approach is not only 
to include non-adiabatic effects in the formation of the 
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FIG. 2: (Color online) Isodensity surfaces at half the satu¬ 
ration denity y9o/2 = 0.08 fm“^ as a function of time for the 
three modes : scf (left), aef (middle) and sef (right). The time 
step between two images are At = 0.675 zs, 1.8 zs and 1.08 zs 
for the scf, aef and sef modes, respectively. 


fragments, but also to provide a well defined value of the 
TKE. Here, the TKE is computed from Coulomb and 
kinetic energies in post scission configurations following 
the method described in Ref. [18]. The TKE are found 
to be 238 MeV, 185 MeV and 163 MeV for the scf, aef, 
and sef modes, respectively. We checked that these re¬ 
sults do not depend on the initial configuration in a given 
valley as long as the initial configuration is not too close 
to the scission point. These results are compared with 
experimental data [26] in Eig. The compact sym¬ 
metric mode is located near the main peak of the TKE 
distribution. This observation is in agreement with the 
experimental data which have attributed this high TKE 
peak to symmetric fission [26]. The lower TKE tail is 
mostly attributed to the asymmetric mode. Our calcula¬ 
tions also predict that the sef mode leads to a low TKE 
where only few events have been observed experimentally. 
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FIG. 3: Experimental distribution of TKE from Ref. [2^. 
The arrows correspond to the mean value of the TKE from 
the three TDHF+BGS calculations. The fragment mass dis¬ 
tribution is shown in the inset. 

A particularly interesting feature of quantum micro¬ 


scopic approaches is the possibility to investigate the role 
of shell effects in the dynamics. Eor instance, TDHE cal¬ 
culations have recently shown the importance of shell ef¬ 
fects in the formation of fragments in heavy-ion collisions 
[SI [35]. Here, shell effects in the tin region (due to the 
proton magic number Z = 50) are expected to be present 
in the symmetric fission mode of ^^^Em [26] |36]. This is 
compatible with the spherical shape of the fragments in 
the scf mode (see Eig. as well as with the high TKE 
associated to this mode. Indeed, magic fragments are dif¬ 
ficult to excite and deform and, thus, fission occurs faster 
as less dissipation is involved, leading to a larger TKE. 
This is also compatible with the short time associated 
with the non-adiabatic descent of the potential to fission 
for the scf mode (see Eig. [^. 

Another possible signature of shell effects in ^^^Em fis¬ 
sion is the narrow peak in the fragment mass distribu¬ 
tion at symmetry (see inset in Eig. [28]. In order to 
see the influence of shell effects on the distributions, we 
have computed the proton number Z and neutron num¬ 
ber N distributions in the fragments at the end of the 
TDHE-fBCS calculations using particle number projec¬ 
tion techniques 1211EZ] with the pfafflan calculated using 
optimized algorithm [38] . The resulting distributions are 
shown in Eig. [^ We clearly see that the Z distribution is 
much sharper for the scf peak than for the other peaks, 
in good agreement with the expectation that this mode 
is dominated by spherical shell effects at Z = 50. The 
peaks in N distributions are all of similar widths, indi¬ 
cating that no shell effect are contributing for neutrons. 
Interestingly enough, we observe a strong odd-even ef¬ 
fect in the N distribution of the scf mode due to neutron 
pairing correlations which is lower in the other modes. 

It should be noted that the Z and N distributions 
shown in Eig. [^ are not expected to reproduce the widths 
of the experimental distribution for two reasons. The first 
reason is that these calculations account only for fluctu¬ 
ations acquired during the non-adiabatic phase. Eluctu- 
ations in the adiabatic phase are indeed expected to be 
important as shown, e.g., in calculations based on the 
time-dependent generator coordinate method [7] or us¬ 
ing the Langevin equation [39]. The second reason is 
that the TDHE+BCS approach is expected to under¬ 
estimate fluctuations of one-body observables @01 HI]. 
Yet, it would be interesting to have a better estimate 
of the fluctuations acquired in the non-adiabatic phase. 
This can be achieved using beyond mean-field theories 
such as stochastic approaches @2] or the time-dependent 
random-phase approximation (TDRPA) @3]. Realistic 
applications of the latter to nuclear dynamics have been 
recently achieved without pairing interaction [44f[46] . 
Solving numerically the TDRPA equations with pairing 
correlations is beyond the scope of this work. Neverthe¬ 
less, we have estimated the enhancement of fluctuations 
due to beyond-mean-fleld effects in the non-adiabatic 
phase of fission for the ^^^Em nucleus fissioning into 
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FIG. 4: Proton (a) and neutron (b) number distributions in 
the fragments for the scf mode (blue dashed line), the aef 
mode (green doted line) and the sef mode (solid red line). 

two doubly magic nuclei in which pairing correla¬ 

tions can be neglected. Numerical details for solving the 
TDRPA equations can be found in [44H47] . As a result, 
the standard deviation for the distribution of the total 
number of nucleons A in the fragment is (Ttdrpa = 2.35 
for the scf mode, to be compared with the TDHF re¬ 
sult (Jtdhf = 1-35. This clearly indicates that beyond 
mean-field fluctuations also play an important role in the 
non-adiabatic phase of fission. 

A fully microscopic approach to the fission process has 
been presented. The path to fission is divided into a slow, 
adiabatic evolution across the fission barrier, followed by 
a faster, non-adiabatic descent of the potential down to 
scission described with a time-dependent mean-field ap¬ 
proach. The method includes pairing correlations and 
can then be applied to superfluid systems across the nu¬ 
clear chart. Application to the fission of ^^^Fm shows 
a good agreement with experimental data. In particu¬ 
lar, this approach can be used to determine the total 
kinetic energy of the fragments without making any as¬ 
sumption on the scission configuration. Quantum shell 
effects are shown to play an important role in the dynam¬ 
ics, in particular in the formation of the fragments in the 
non-adiabatic phase. The method could be generalized 


to systems with finite temperature in order to study the 
disappearance of shell effects [48l [49] . Beyond mean-field 
fluctuations are shown to be important. Thus, a quanti¬ 
tative description of fragment distributions requires fur¬ 
ther developments of beyond mean-field approaches. 
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